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ABSTRACT 

The plasma /3 (the ratio of the plasma pressure to the magnetic pressure) of a system can have a 
large effect on its dynamics as high /3 enhances the effects of pressure anisotropics. We investigate the 
effects of j3 in a system of stacked current sheets that break up into magnetic islands due to magnetic 
reconnection, which is analogous to the compressed heliospheric current sheet in the heliosheath. We 
find significant differences between systems with low and high initial values of /?. At low /? growing 
magnetic islands are modestly elongated and become round as contraction releases magnetic stress and 
reduces magnetic energy. At high f3 the increase of the parallel pressure in contracting islands causes 
saturation of modestly elongated islands as island cores approach the marginal firehose condition. Only 
highly elongated islands reach finite size. The anisotropy within these islands prevents full contraction, 
leading to a final state of highly elongated islands in which further reconnection is suppressed. The 
elongation of islands at finite j5 is further enhanced by reducing the electron-to-ion mass ratio, to 
more realistic values. The results are directly relevant to reconnection in the sectored region of the 
heliosheath where there is evidence that elongated islands are present, and possibly to other high /? 
systems such as astrophysical accretion flows and the magnetosphere of Saturn. 



1. INTRODUCTION 

At the outer edges of the solar system, the solar wind 
pushes up against the interstellar medium. Since the so- 
lar wind is moving at supersonic speeds, it forms a shock 
known as the termination shock. The region between 
the termination shock and the interstellar medium is re- 
ferred to as the heliosheath. Flapping of the heliosp heric 
curre nt sheet produces sectored magnetic fields ([Smith! 
2001) that have a significant latitudinal extent. There 
has been research suggesting that the current sheets be- 
tween the sectored fields found in the heliosheath are 
compressed to the point that collisionless reconnection 
begins t o occur, resulting i n the formation o f magnetic 
islands ([Drake et al.l 120101 : iCzechowski etaLl feOlQ). A 
turbulent magnetohydrodynamic (MHD) model of the 
reconn ection of the sectored fields has also been pro- 
posed ([Lazarian fc Opherl l2QQ9h although we will argue 
later that the Voyager data are inconsistent with this 
hypothesis. 

An important question, however, is whether the 
conventional trea tment of collisionless reconnection 
(Sha v et al.l [2007) is valid in the heliosheath, where it 
was suggested that the pick-up ion (PUI) population in- 
crea ses the plasma pressure compared with values at 1 
AU (|ZankHl99arRichardson etaLll2QQa IWu et al.ll2009h . 
Although both Voyager spacecraft are currently taking 
data in the heliosheath, the energy range of the de- 
tectors does not cover the PUIs, so it is difficult to 
make a reliable estimate of the value for /3, the ra- 
tio of the plasma pres sure to the magnetic pressure 
([Richardson et al.ll2008[ ). Global MHD simulations sug- 
gest, however, that j3 varies from 8 to 0.5 between the 
termination shock and the interstellar medium with the 
highest /3 just do wnstream of the termination shock 
([Drake et al.ll2010] ). Although this simulation does not 
include a separate pick-up ion population, it provides a 



rough estimate for the expected values for f3 and mo- 
tivates the range of j3 in our study. In this study we 
investigate the impact of j3 on the dynamics of reconnec- 
tion and the formation of magnetic islands relevant to 
the sectored heliosheath. 

To begin reconnection and island formation a current 
sheet needs to be compressed to approximately the ion 
inertial scale d{ = c/co> p i, where c is the speed of light 
and u; D i is the ion plasma frequency, ([Cassak et al. 2005; 
lYama da 2007). At this point the current sheet becomes 
unstable to the collisionless tearing mode. Upstream of 
the termination shock the heliosp heric curren t sheet has 
a thickness of around 10, 000 km (Smith 2001) and is pre- 
dicted to compress to around 2500 km just downstream 
of the shock. In the upstream region the plasma den- 
sity measured at Voyager 2 is around 0.001cm -3 , cor- 
responding to an ion skin depth of around 7, 200 km. 
In the downstream region the density is compressed to 
around 0.003cm -3 , corresponding to an ion skin depth 
of around 4200 km. Thus, the compression of the cur- 
rent sheets across the termination shock should trigger 
collisionless reconnection in the heliosheath and some 
Voyager 1 and 2 observations support this hypothesis 
(jOpher et al.ll2Qlih . 

In our system we will be examining symmetric current 
sheets with no guide field (init ial out-of-plane m a gnetic 
field). In an analytic study iBrittnacher et al.l (p.995) 
showed that when pi/wo w 1, the fastest growing linear 
mode occurs at kwo ~ 0.5, where k is the wavenumber 
of the tearing mode, w$ is the half width of the current 
sheet, and pi is the ion gyroradius. Since pi = \ffTidi, 
where fy is the plasma beta based on the ion pressure, 
the current sheet thickness is comparable to that in our 
simulation. 

Due to the collisionless nature of the plasma, pres- 
sure anisotropics (P\\ ^ P_l, where || and _L are de- 



termined with respect to the magnetic field) can form 
during r econnection. Ferm i acceleration in contracting 
islands (jDrake et al.l [2006) and adiabatic cooling as B 
decreases due to conservation of the magnetic moment 
/i oc v\jB oc P±/B both drive P\\ > P±. When an 
anisotropy is formed with Py > P± the tension in bent 
magnetic fields weakens. The fluid momentum equation 
with an anisotropy becomes 
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where p is the mass density of the plasma, v is the bulk 
velocity, and B is the magnetic field. For j3\\ = /3± the 
pressure equation reduces to the standard MHD equa- 
tion. When /3\\ > f3± the tension force is reduced. At 
high /3 this reduction of tension force is noticeable for 
even slight anisotropics in the pressure. For f3\\ — f3± 
large enough, the tension force drops to zero, or even 
becomes negative. Since the tension of field lines acts 
as a restoring force for Alfven waves in standard MHD, 
the negative sign causes this oscillation to become an in- 
stability known as the firehose instability (Parker 1958) 
for: 

£|| - P± > 2- (2) 

This instability is fueled by the free energy contained in 
the pressure anisotropy. The firehose instability causes 
magnetic field lines to kink, which eventually relieves the 
pressure anisotropy by causing scattering. 

Alternatively, when f3\\ — f3± is negative and large 
enough in magnitude, other instabilities can occur. The 
mirror-mode instability and the ion cyclotron instability 
both occur when P± > Pu. For larger /3m the mirror- 
mode becomes unstable at smaller values of \/3\\ — /3±\ 
than the ion cyclotron mode, so the marginal mirror- 
mode criterion acts as the boundary between the stable 
and unst able regions. Ba sed on fluid theory assuming 
T e = Ti (Hasegawa 1969), the mirror- mode instability 
occurs when 



P± ~ P\\ > 



P±' 



(3) 



There are also kinetic modifications that can be made 
to the marginal instability criteria for firehose, mirror- 
mode, and ion cyclotron which make them more accu- 
rate. Although a rigorous analytic theory is not avail- 
able, ther e are models that approximate the in stability 
very well ([Hellinger et al.l 120061 : IBale et all 12009( 1 . How- 
ever, for simplicity we will just consider the conditions 
based on fluid theory. 

In this study, we simulate several stacked current 
sheets similar to the compressed sectored heliospheric 
fields and associated current sheets, and follow the devel- 
opment of reconnection and islands. We implement this 
system in a two-dimensional particle-in-cell (PIC) code 
and vary the temperature of the background plasma to 
test the dependence on f3. We observe that in finite j3 e 
systems (/3 e > 0.5), very elongated islands form as op- 
posed to the modest-aspect-ratio islands found at low j3 e 
(/? e < 0.5), where f3 e is the /3 based on the electron pres- 
sure. At high /?, the increased P\\ due to the Fermi re- 
flection of electrons within islands saturates the normal 
modest-aspect-ratio islands. Fermi reflection in highly 



elongated islands is less efficient because of the increased 
bounce time of the electrons so these islands are able to 
reach finite amplitude. At late time, however, even these 
elongated islands exhibit anisotropy instabilities, from 
Fermi reflection of both ions and electrons. As a result, 
late-time magnetic islands remain highly elongated and 
do not become round as in the low-/3 regime. This result 
has significant implications for the structure of islands 
that would be measured in the heliosheath. Although /3 e 
is rather moderate in the heliosheath, we find a mass ra- 
tio dependence suggesting long islands for a broad range 
of j3 e in realistic mass ratios. A large (3 however may be 
necessary to sustain the elongation of these islands. 

2. COMPUTATIONAL MODEL 

Our simulations are performed with the PIC code p3d. 
The in itial conditio ns consist of eight Harris current 
sheets ([Harrisl |1962[ ) where the magnetic pressure bal- 
ances the plasma pressure. Each Harris sheet consists 
of a taiih.(y/wo) and sech (y/wo) profile along the y di- 
rection for an x-directed magnetic field and the density, 
respectively. The peak density of the Harris sheets is no- 
In addition, there is a uniform background population 
that has a density of n^ = 0.2no- These simulations are 
done in two dimensions so d/dz = 0, where z is out-of- 
plane, parallel to the initial current. 

The code uses normalized units. The timescale is 
normalized to the ion cyclotron time ft^ 1 . The dis- 
tance scales are normalized to the ion inertial length 
di = c/cjpi. Thus, the velocity is normalized to the 
Alfven speed va- The magnetic field is normalized to 
the asymptotic value of the reversed magnetic field Bq. 
The density is normalized to no- The pressure is nor- 
malized to Po = n§miv\ — Bq/4tt. The temperature is 
normalized to T = rriiv\. 

In order to vary the /3 of these simulations we vary 
the temperature of the background population TJ,. This 
background temperature is the same for both ions and 
electrons. The Harris equilibrium is used to balance 
the sharp change in the magnetic field strength across 
the current sheets, while the background represents the 
PUIs and has the greatest influence on late-time recon- 
nection dynamics. We performed simulations for /? = 0.2, 
1, 2, 3, and 4.8, where /3 is based on the pressure 
in the asymptotic field with density n^. Each simula- 
tion was advanced for a time of 120^ with a time 
resolution dt = 0.0040^ . The simulations are on a 
204. Sdi x 102. Adi domain with a grid scale resolution of 
A x = A y = 0.05di. In order to complete such large runs, 
unless otherwise specified, we used 25 for the mass ratio 
of the ions to electrons. This makes it easier to resolve 
small electron scales. In order to lessen the separation 
between the field and particle timescales, we set the ratio 
of the speed of light to the Alfven speed, c/ca, to 25 (in 
the heliosheath a more realistic value is near 6000). Re- 
connection is insensitive to the value of c/ca- We start 
with a half-thickness for the current sheet wq = 0.5<ii, 
so that collisionless reconnection can begin from particle 
noise. The temperature in the Harris sheet is 0.25Xb for 
both ions and electrons, and there is no guide field. The 
largest j3 we simulated was 4.8 since the electron thermal 
velocity v t he ~ 0.7c. Larger j3 would begin to have sig- 
nificant unphysical relativist ic effects, due to our lowered 
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Fig. 1. — Normal component of the magnetic field, B y , for (3 = 2 
at t = 1517^ . The bipolar signatures in the current sheet indicate 
the presence of x-lines arising from the collisionless tearing mode. 

ratio of c/ca- 

The simulation does not precisely describe the he- 
liosheath but illustrates important physics that should be 
found there. The ion pressure in the heliosheath is much 
larger than the electron pressure, and thus (3 e is actually 
quite moderate. However, we will show that at realistic 
mass ratios we still expect elongated islands. The large 
/?, mostly due to ion pressure, is however significant be- 
cause it allows the elongation to persist. We do not have 
a separate population of PUIs, and the magnetic field 
configuration is a Harris sheet rather than the rotated 
fields (where \B\ is constant in a cut through the current 
sheet) found in the heliosheath. Future simulations tai- 
lored to the specific parameters of the heliosheath may 
help our understanding of this phenomenon. 

3. RESULTS 

The early development of a run with f3 = 2 is shown in 
FigureHJ Not surprisingly a wave mode with k x wo ~ 0.5 
clearly emerges. The finite f3 background plasma does 
not have a strong effect on the wavelength of linear tear- 
ing. During this time, within the current sheets, an 
anisotropy in the electron pressure begins to develop with 
P e \\ > P e ±. The electrons moving at the thermal velocity 
are able to bounce between the two ends of the islands 
which have lengths of ~ 6^. Comparisons with runs at 
different f3 show that at this early time, the lengths of 
the islands appear to be insensitive to f3. 

There are two important timescales controlling the dy- 
namics: the time it takes for the ions to accelerate to 
Alfvenic outflow speeds and the time it takes for a signifi- 
cant electron pressure anisotropy to develop. If electrons 
bounce several times between the two ends of a contract- 
ing island, an anisotropy develops which approaches the 
firehose instability boundary. This is because the bounc- 
ing electrons gain energy in the parallel direction. The 
time for an electron to bounce off the edge of an island 
and then return to its original position is thus a measure 
of the time for significant anisotropy to develop. 

The tearing instability is driven by the tension in the 
newly reconnected magnetic fields. Since anisotropics 
cause a weakening of the magnetic tension, the tearing 
mode can be suppressed by strong anisotropy within an 
island. 

Reconnecting magnetic field lines, by relaxing their 
tension, accelerate ions up to Alfvenic speeds. If several 
bounces occur during the time required for ions to be 
accelerated up to the Alfvenic outflow speed from the x- 
line, the developing anisotropy slows the ion outflow and 
essentially stops the growth of the tearing mode. How- 
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Fig. 2. — Ion outflow velocity, vi x , vs. the position x, for t = 
2bQ^ 1 (green), 30^ 1 (red), and 3bQ~ { 1 (black) for the (3 = 2 run. 
To reduce noise we do a five-point smoothing of vi x in both the 
x and y directions. The blue curve is a line of slope 0.1, which 
corresponds to a convective growth time t a of around 10f2^ . 

ever, since the bounce time is proportional to the length 
of the islands, the growth of sufficiently long wavelength 
tearing modes can continue. 

Near an x-line adjacent to a growing island, the outflow 
velocity of ions, to first approximation, linearly increases 
with distance (see Figure[5]) as 



1 

Vix — —X. 



(4) 



At this point t a is just defined as the inverse slope of 
the relationship between v\ x and x. Since according to 
Figure[2]the slope is constant, this implies that ions ac- 
celerate away from the x-line in an exponential fashion. 
By integrating, the time, t, for an ion to accelerate from 
the initial position, #o, to the final position, Xf, can be 
obtained. 
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= rr^t/ta 



Xf = x$e 



(5) 



(6) 



An approximate measure for the characteristic timescale 
for acceleration away from the x-line up to the Alfven 
speed is t « t a: the acceleration time. As seen in Fig- 
ure[2l the acceleration time at t — 30^ 1 is of order 
~ lOfi^ . This acceleration time is approximate, can 
vary by a factor of as much as two, and appears to be 
insensitive to /?. 

The bounce time can be estimated based on the ther- 
mal velocity of the electrons, v t he, and the length of the 
island, L: 

t b = — = -\[^, (7) 

^the v A V Pe rrii 

where /? e is the /? determined solely from the plasma 
pressure derived from the electrons. 




Fig. 3. — Regions where the plasma anisotropy causes instabil- 
ity. White regions are unstable to the firehose, black regions are 
unstable to the mirror-mode, and red are stable. The green lines 
are magnetic field lines. This plot shows one current layer taken 



from the (3 = 2 run at (a) t = 25Q~ { x and (b) t = 40Q~ 
aspect ratio is distorted to make the islands more visible. 
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Equating the empirical acceleration time t a = lOQ^ 1 , 
and the bounce time £5, Equation ([7]) , a critical island 
length can be found: 



Lcrit ~ lOdi* p t 






(8) 



For islands with L < L cr n, the anisotropy will stop the 
tearing instability. Islands smaller than L crit can still 

form, but t hey quickly satura te. A similar saturation 

was found in lKarimabadi et all (j2QQ5h . However, in their 
simulations the size of the computational domain was 
12.6pi and L cr n = lOOp^, where pi is the ion Larmor 
radius. Thus, the development of long wavelength islands 
was not observed. 

For the case of /3 = 2 (/3 e = 1) and L « 6di, 
tb ~ i.2^ 1 . This time is much less than the accel- 
eration time, so there is enough time for a significant 
anisotropy to develop before a significant x-line is estab- 
lished. This anisotropy can be seen in Figure^a), which 
shows the regions from the /3 — 2 run that are unstable 
to the firehose instability. The unstable regions occur 
inside the islands and stop further growth of the short 
wavelength tearing modes. The islands that continue to 
grow correspond to longer wavelength, with L w 40di 
and tb ~ 80^ « t a . Thus, the anisotropy develops 
slowly enough for reconnect ion to develop. This can be 
seen in Figure[3fb). 

As can be seen in FigureHJ by £ = 510"^, /3 has a 
significant influence on the structure of islands. The is- 
lands for /3 = 0.2 have much shorter wavelength than for 
/3 — 2 and 4.8. In other words, there are more locations 
where reconnection proceeds in the case of low j3. This 
phenomenon is expected based on the previous analysis, 
Lcrit cx \fWe- Since L cr i t is proportional to the square 
root of the mass ratio yjmijm e , we expect to find much 
longer islands in the real mass ratio limit. To test this, we 
perform a fi — 0.2 simulation with mi/m e = 100. In this 
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Fig. 4. — Out-of-plane current, j z , at t = 51^ ci 1 for /3 of (a) 0.2, 
(b) 2, and (c) 4.8. 




Fig. 5. — Out-of-plane current, j z , for (3 = 0.2 at (a) t = 40^ ci 1 

for mi/m e = 25 and (b) t = 60f2^T for mi/m e = 100. The aspect 
ratio is distorted to make the islands more visible. 

case, we reduce the ^-domain by a factor of four with re- 
spect to FigurelUa), examining only two current sheets. 
We double the resolution in order to resolve the small 
electron scales and reduce the ratio of the speed of light 
to the Alfven speed to 15. There is a clear dependence on 
rrii/m e shown in Figure[5j where we compare the bottom 
two current sheets of Figure|2^a) to the new simulation. 
We find the islands to be significantly longer, confirming 
our prediction. Since mi/m e ^> 100 in the heliosphere, 
long islands are almost always expected, unless f3 e is very 
small. 

As elongated islands grow at high /3, anisotropics 
within them also develop even though the anisotropics 
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Fig. 6. — Ion density, ni, for (3 
between t = 41 and 66f2~T . 
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Fig. 7. — Fraction of grid points unstable to the firehose insta- 
bility vs. time for (3 = 0.2, 1, 2, 3, and 4.8. 
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Fig. 8. — Two-dimensional histogram of anisotropy (a = P±/P\\) 
vs. (3\\ for (3 = 2 at times from top to bottom t = 80, 120, 

and 160^2^ . The blue line represents the marginal condition for 
mirror-mode instability. Points above this curve are unstable. The 
red line represents the firehose marginal stability condition. Points 
below this curve are unstable. The color bar represents the number 
of points with a particular a and j3u . 



do not suppress island growth. The anisotropy sur- 
passes the firehose condition in the center of the islands. 
These anisotropics are lik ely caused by the Fermi mech- 
anism (Drak eet al.l 12006). The dynamics of this accel- 
eration mechanism will be discussed in a future paper. 
The contraction of islands can be seen in Figure[6l The 
higher density regions inside of the islands move inward 
at Alfvenic speeds. At around t = 61 — 66^"^ the islands 
begin to kink, which indicates the onset of an anisotropy 
instability. 

The short-wavelength mode is caused by the tem- 
perature anisotropy due to the outflow from the x- 
line streaming through the plasma entering the exhaust 
across the separatrix and the Fermi acceleration of elec- 
trons bouncing in the island. Based on the similarities in 
growth rate and other signatures that will be discussed 
in a future paper, this mode appears to be associated 
with the Weibel instability. 

The anisotropics that develop during the reconnec- 
tion simulation do not grow without bound. In Fig- 
ure[7J we plot the fraction of grid points that are un- 
stable to the firehose instability. As time advances and 
the anisotropics begin to form, the number of grid points 
unstable to the firehose instability increases. However, at 
t ~ 800^, the number of unstable grid points begins to 
saturate. Since it takes place soon after the onset of the 
kinking of the islands, the saturation is likely because 



the anisotropy is reduced via scattering by the Weibel 
and firehose instabilities. Additionally the saturation oc- 
curs soon after the unreconnected flux is exhausted. By 
60^ 1 the islands have grown enough so that the islands 
on adjacent current sheets begin to interact. This is an 
additional reason for the saturation of the firehose unsta- 
ble area: there is no more space into which the firehose 
unstable islands can expand. 

At late time, the anisotropy of the system is confined 
within the boundaries of the marginal firehose (Equa- 
tion©) and mirror- mode instabilities (Equation©) in 
a man ner similar to that seen in observations of the solar 
wind (jHellinger et al.l 120061 : [Bale et al.ll2QQ9h and in ear - 



lier low- /3 current sheet simulations (jDrake et al.l 12010). 
Figure[8] shows the data for our system in the space of 
(a:,/3||) where a = P±/P\\. This plot is generated by cal- 
culating the anisotropy a and the /3y for each grid point. 
The plot is a two-dimensional histogram of grid points 
in (a,/?n) space, where /3m is calculated based on Py. The 
parallel and perpendicular pressures are calculated by 
taking the diagonal components of the pressure tensor 
after rotating into the frame of the local magnetic field, 
such that the two perpendicular components are equal. 
We look at the distribution at t = 80, 120, and 160^ . 
At early times the anisotropics have not yet fully de- 
veloped and the plasma still occupies a small region in 
(ck,/3||) space. By t = 1201^ the anisotropy has reached 
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Fig. 10. Distribution of A at t = HOfi^ 1 for the (a) = 0.2 
case and (b) /3 = 4.8 case. The dotted lines are at A = 90° and 
270° where we expect to find peaks in the distribution. 



Fig. 9. — Out-of-plane current, jz(a), and stability(b) at t = 
120r2^ 1 for the (3 = 4.8 case. White regions in part (b) are unstable 
to the firehose, and red are stable. 

the two stability boundaries, and continues to be con- 
fined between these two boundaries at t = 160^"/, even 
as the average j3 increases. The anisotropy reaches the 
stability boundaries at a time after the short wavelength 
Weibel modes have dissipated. Since at this point there 
are no longer large regions with essentially zero mag- 
netic fields, the firehose and mirror- mode instabilities 
are what determine the boundaries of the temperature 
anisotropics. There are no clear signatures of the clas- 
sical mirror- mode instability at this time. The firehose 
and mirror-mode instabilities may be hard to distinguish 
among the turbulent interacting magnetic islands, or the 
islands may just stop generating anisotropy as they ap- 
proach the instability boundaries. 

The islands maintain an elongated form for the simu- 
lation shown in Figure|4^c) clear until t = 120^"^, the 
latest time simulated for f3 = 4.8. This is shown in Fig- 
ure[9fa) showing the out-of-plane current for t = 120^"^. 
Since the edges of the islands are pushing against the fire- 
hose instability, the tension force in the magnetic fields is 
eliminated. This can be seen in Figure[9jb) which shows 
the regions that are unstable to the firehose instability. 

4. CONCLUSIONS 

The magnetic islands that reach a significant ampli- 
tude are much more elongated at high /3 e than at low j3 e . 
These elongated islands should be found even for moder- 
ate values of /3 e at realistic mass ratios. Island elongation 
is caused by the suppression of the shorter wavelength 
tearing modes by pressure anisotropics (Py > P±) that 
develop due to the Fermi acceleration of electrons. Later 
in time the plasma develops pressure anisotropics of both 
ions and electrons that are limited by the firehose and 
Weibel instabilities. A Weibel mode develops that kinks 
the magnetic field lines. In the regime with a real mass 
ratio we would expect even longer islands to form, where 
multiple wavelengths of the firehose instability could de- 
velop. At late time the fraction of points unstable to 
the firehose instability saturates, and the anisotropy is 
confined between the mirror-mode and firehose instabil- 



ity boundaries. The long islands persist due to the low 
requirement of anisotropy to reach the marginal firehose 
condition at high j3. For even small anisotropics the ten- 
sion in the magnetic fields is removed. 

When encountering magnetic islands in the he- 
liosheath, we predict the formation of similar extended, 
sausage-shaped islands rather than the more round is - 
lands found in low-/3 simulations (Drak e~et al.l [2010). 
The cores of these islands should also be at the marginal 
firehose condition, so the magnetic tension that drives 
them to become round vanishes. We would thus expect 
these sausage shapes to persist long after the islands have 
ceased growing, and thus could be found even in regions 
where reconnection is no longer occurring. 

These elongated islands exhibit signatures that can be 
seen in Voyager data. In particular, Voyager measures 
all three components of the magnetic field. Of particular 
interest for the explorations of islands that grow in the 
ecliptic plane is the angle A = tan -1 (Bt/Br), where Bt 
and Br are the azimuthal and radial magnetic fields, re- 
spectively. A = 90° and 270° correspond to the azimuthal 
unreconnected sectored heliosheath magnetic fields. De- 
viation of A from 90° and 270° indicates some process 
is distorting the sectored field. Voyager data show the 
distribution of A is peaked in the two azimut hal direc- 
tions, A = 90° and 270° (jQpher et all 1201 If ). These 
peaks are significantly broader in the heliosheath than 
upstream, indicating that reconnection or another mech- 
anism is disturbing the heliosheath field. The observed 
Voyager distributi on of A is consistent with that found in 
high-/3 simulations (jQpher et al.ll20Tll ). Since the islands 
are elongated, the magnetic fields tend to remain primar- 
ily in the azimuthal direction even well after the islands 
begin to interact with each other. Round islands, such 
as would be expected from an MHD model or a low j3 
kinetic model, are not consistent with observations since 
they produce much broader A distributions. Thus, MHD 
reconnection (Lazarian & Opher 2009) in the heliosheath 
seems to be ruled out. Shown in FigureflOlis the distribu- 
tion of A from the simulations at j3 = 0.2 (Figure[H(a)), 
and ft = 4.8 (Figured^)) at t = HOfi^ 1 . The high /? 
simulation which has elongated islands retains the two 
peaks at A = 90° and A = 270°. The long islands have 
a larger magnetic field in the azimuthal direction than 



the radial, resulting in peaks in the A distribution, but 
the shorter islands become round having a magnetic field 
with similar strength in both directions, resulting in a 
broad distribution in A. The loss of tension in a finite /3 
plasma prevents the complete release of magnetic energy 
that would be expected in an MHD model. A complete 
understanding of the j3 dependence of magnetic islands 
is essential in order to obtain reliable signatures that can 
be compared with Voyager data. 

In this work there was no out-of-plane guide magnetic 
field. In the heliospheric current sheet, the magnetic field 
rotates from one direction to the other keeping a con- 
stant magnitude rather than passing through zero ([Smith! 
|2QQ1[ ) . A guide field would cause the center of the islands 
to have a much lower (3 since the magnetic field does not 
go to zero. Because of this magnetic field, we would not 
expect the Weibel instability to develop. In real systems 
there is frequently a guide field, so this would be worth 
further investigation. 

In astrophysical accretion disks, reconnect ion plays a 
role in determining the s aturation of the magnetorota- 
tional instability (MRI) dSano et al.ll2QQ4f ). The satu- 
ration of MRI is strongly dependent on the dissipation 
of the magnetic field due to reconnection. For the high 
j3 in accretion disks, suppression of the most strongly 
growing small islands may significantly impact the sat- 
uration of the MRI. Since j3 is typically larger than 100 
in these structures, the only surviving islands would be 
so long that it is likely that much of the magnetic free 
energy would not be dissipated. Further, since the MRI 



requires magnetic tension, the absence of tens ion could 
limit the development of the instability. iSharma et al.l 
(2006) perform a simulation showing an enhancement of 
the growth of MRI due to anisotropics with P± > P\\ , 
which enhances the magnetic tension, caused by \i con- 
servation as a magnetic field develops. They do not cap- 
ture the physics of reconnection and Fermi acceleration 
in magnetic islands that would generate anisotropics with 
P± < P||, which removes magnetic tension. These two 
competing sources of anisotropy, both affect the tension 
and thus the growth of the MRI. The relative importance 
of these mechanisms needs to be explored. 

Reconnection at high /3, although relatively rare in the 
terrestrial magne t osphere, is also foun d in the magneto- 
sphere of Saturn (Ma sters et al.ll2011[ ). Magnetic islands 
were discovered in a region where f3 is larger than 10. 
The /? dependence of the growth of finite-sized magnetic 
islands may lead to a better understanding of these find- 
ings. 

Since the development of elongated islands requires 
only moderate /?, we expect to see the development of 
longer islands than expected in lowe r systems such 
as the magnetosphere. In contrast to Karimabadi et al.l 
(2005), these longer islands can grow to a large enough 
size to play a role in magnetospheric dynamics. Since 
the j3 of the magnetosphere is not exceptionally large it 
is unlikely that the persisting anisotropy is enough to 
keep the islands from eventually becoming round. 

The computations were performed at the National En- 
ergy Research Scientific Computing Center. 
1,1 Top 
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